How does the geodesic rule really work for global symmetry breaking first order phase 

transitions? 
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The chain of events usually understood to lead to the formation of topological defects during phase 
transitions is known as the Kibble mechanism. A central component of the mechanism is the so- 
called "geodesic rule". Although in the Abelian Higgs model the validity of the geodesic rule has 
been questioned recently, it is known to be valid on energetic grounds for a global U(l) symmetry 
breaking transition. However, even for these globally symmetric models no dynamical analisys of 
the rule has been carried to this date, and some points as to how events proceed still remain obscure. 
This paper tries to clarify the dynamics of the geodesic rule in the context of a global U(l) model. 
With an appropriate ansatz for the field modulus we find a family of analytical expressions, phase 
walls, that accounts for both geodesic and nongeodesic configurations. We then show how the latter 
ones are unstable and decay into the former by nucleating pairs of defects. Finnally, we try to give 
a physical perspective of how the geodesic rule might really work in these transitions. 
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I. INTRODUCTION 

According to the current models in particle physics, 
symmetry breaking phase transitions are expected to 
have occurred in the early stages of the evolution of the 
universe. There are two mechanisms by which these tran- 
sitions may have taken place: either by the formation of 
bubbles of the new phase within the old one (i.e., first or- 
der phase transition), or by spinodal decomposition (i.e., 
second order transition). Although we still do not know 
which one did actually take place in each particular case, 
for the electroweak phase transition for instance common 
opinion inclines more towards the first of these two possi- 
bilities. As the first order transition scenario would have 
it, bubbles of the new phase nucleated within the old 
one (the nucleation process being described by instan- 
ton methods as far as the WKB approximation remains 
valid p]]), and subsequently expanded and collided with 
each other until they occupied all of the available vol- 
ume at the time at which the transition was completed. 
In the process of bubble collision though, the possibil- 
ity arises that regions of the old phase become trapped 
within the new one, giving birth to topologically stable 
localized energy concentrations known as topological de- 
fects (for recent reviews see refs. Q]), much in the same 
way in which these structures are known to appear in 
condensed matter phase transitions. 

From a theoretical point of view, topological defects 
will appear whenever a symmetry group G is sponta- 
neously broken to a smaller group H such that the result- 
ing vacuum manifold M — G/H has a non-trivial topol- 
ogy: cosmic strings for instance (vortices in two space 
dimensions) will form whenever the first homotopy group 
of M is non-trivial, i.e., 7Ti(M) ^ 1. To see how this could 



happen in detail, let us consider the Lagrangian 
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for a complex scalar field <f>. Let us assume that V is 
of the type V = j(\$>\ 2 — ?7 2 ) 2 , and that its parameters 
are functions of the temperature such that at high tem- 
peratures $ = is the only minimum of V , while at 
zero temperature all the |$| =7] states correspond to 
different degenerate minima. Then the structure of the 
vacuum manifold will be that of S . iti(S ) ^ 1 however, 
and thus we can form non-contractible loops in the vac- 
uum manifold. That is, the model admits cosmic string 
solutions. 

The way in which these configurations would actually 
arise during a phase transition is via the Kibble mech- 
anism In the context of a first order transition the 
basic ideas behind it are that bubbles are nucleated with 
random phases of the field, and that when two regions 
in which the phase takes different values encounter each 
other the phase should interpolate between this two re- 
gions following a geodesic path in the vacuum manifold 
-the so-called geodesic rule. A possible scenario for vor- 
tex formation in two space dimensions would then look 
like this: three bubbles with respective phases of 0, 2tt/3, 
and 47r/3 collide simultaneously. Then, if we walk from 
the first bubble to the second one the geodesic rule tells 
us that the phase has to grow from to 2ir/3 along our 
path (as opposed to decreasing from to — 4ir/3). In 
the same manner, walking from the second bubble to 
the third one we will see the phase grow from 2tt/3 to 
47r/3, and if we finally walk from the third bubble back 
to the first one again we will see a change from 47r/3 to 
2n. The phase will have thus wounded up by 2n in our 
path, having traversed the whole of the vacuum manifold 
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once along the way. Continuity of the field everywhere 
in the region inside our path demands then that the field 
be zero at some point in this region, namely the vortex 
core. In the limit in which the bubbles extend to infin- 
ity, outwards from the center of collision, removal of this 
vortex would cost us an infinite amount of energy, since 
it would involve unwinding the field configuration over 
an infinite volume. The vortex is thus said to be topo- 
logically stable. In three space dimensions, the resulting 
object would obviously be a string, rather than a vor- 
tex. Clearly though there are other ways in which strings 
could be formed. Collisions of more than three bubbles 
could also lead to string formation, or, for instance, two of 
the bubbles could collide first, the third one hitting them 
only at some later time while the phase is still equilibrat- 
ing within the other two. This event in particular will be 
far more likely than a simultaneous three way (or higher 
order) collision, and it is probably the dominating pro- 
cess by which strings are formed (especially if nucleation 
probabilities are as low as required for WKB methods to 
be valid) . When two bubbles collide, the phase will try to 
reach a homogeneous distribution within the single true 
vacuum cavity formed after the collision. However, in the 
absence of coupling between the scalar and the plasma 
that would also be present in a cosmological setting, this 
process is never completed -essentially because the ve- 
locity at which the phase propagates inside the bubbles 
is the same as that with which the bubble walls expand. 
Thus, this does not substantially modify the picture of 
defect formation described above. 

The geodesic rule is therefore central in that not only 
gives us a heuristic view of how bubbles interact after 
they collide, but it also yields an estimate of the density 
of defects produced in a phase transition, i.e., about one 
per every three bubbles (times the probability that the 
three bubbles have the correct phases of course). If one 
is concerned about the density of defects that are present 
after a cosmological phase transition, the question of the 
validity of the rule is therefore of obvious importance. 
Despite its importance however, the rule still presents 
serious problems in its understanding. For globally sym- 
metric models the role of the geodesic rule was studied 
by Srivastava some time ago. In M he established its 
validity within these type of models mainly on energetic 
grounds, the geodesic path being the least energetic of 
all the possible paths that the phase can follow. Re- 
cently however a number of authors have argued for and 
against the validity of the rule § , or observed what seems 
to be a breakdown of the rule in numerical simulations 
H , within the context of the Abelian Higgs model. In 
this case it is rather less clear whether one can reason 
on energetic grounds since the phase difference between 
two points is not a gauge invariant quantity by itself. In 
principle therefore one could always choose a gauge in 
which the phase is zero everywhere, and it would not be 
immediately obvious how the geodesic rule would work 
in this situation. The goal of this paper is to take a step 
back and try to give a dynamical account of how the 



geodesic rule may actually work for a globally U(l) sym- 
metric model. In order to do this we will use a type of 
analytical solutions, hereafter referred to as phase walls, 
by which the phase interpolates between its values within 
two bubbles. We will also resort to numerical simulations 
to study the stability of different phase wall configura- 
tions, and finally we will discuss the physical settings in 
which non-geodesic phase walls can appear. The paper is 
organized as follows: in section II we give a brief account 
of the model and the results previously found in Q rel- 
evant for this paper. We then extend our treatment in 
section III to include the case in which the scalar field is 
not dissipatively coupled and we have 2+1 dimensions. 
In section IV we study the stability of non-geodesic phase 
walls. In V we study the influence of our work in three 
bubble collisions and topological defect nucleation, and 
in VI we present a general view of how the geodesic rule 
might work in globally symmetric models. Finally, the 
conclusions and the implications that our analysis may 
have for the Abelian Higgs mode are presented in section 
VII. 



II. 1 DIMENSIONAL PHASE WALLS IN A 
DAMPING ENVIRONMENT 

A. expression for the interpolating phase 

Consider the Lagrangian (Q) for a complex field <f>. We 
will use the same form of potential that was used in 
that is, 



V = A 



^(l*|-^) 2 -^l*| 3 



(2) 



This is just a quartic potential with a minimum at |$| = 
(the false vacuum), and a set of minima connected by 
a U(l) transformation (true vacuum) at |$| = p tv = 

^(3 + e + \J (3 + e) 2 — 8), towards which the false vac- 
uum will decay via bubble nucleation. It is the dimen- 
sionless parameter e that is responsible for breaking the 
degeneracy between the true and the false vacua. 
The equations of motion for this system are then 



<9 M ^$ = -dV/d§. 



(3) 



For the potential (g), approximate solutions of (g) exist 
for small values of e, the so-called thin wall regime 
and are of the form 



1*1 
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l-tanh(^( X -iio) 



(4) 



where Rq is the bubble radius at nucleation time and 
X 2 = | x | 2 — t 2 . The bubble then grows with increasingly 
fast speed and its walls quickly reach velocities of order 
1. In M however the authors were concerned about a 
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model with overdamped motion of the walls due to the 
interaction with a surrounding plasma. In order to model 
this effect, they inserted a frictional term for the modulus 
of the field in the equation of motion, namely 



|i 0"* + 7|*| 



-dV/d&, 



(5) 



where |$|= d\$\/dt,6 is the phase of the field, and 7 
stands for the friction coefficient -which will serve as pa- 
rameter under which we will hide our lack of knowledge 
about the detailed interaction between the wall and the 
plasma. It can be shown (see § ) that this equation does 
indeed posses a solution that shows the desired type of 
overdamped motion for the bubble walls. In the thin wall 
limit, this solution can be written as 
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1 — tanh 



V\rj (r - v ter t - R ) 

t 2 yrr^r; 



(6) 



which is simply a spherically symmetric Lorentz- 
contracted moving domain wall with a velocity w ter of 
the form Vter ~ /jPt v , where 5 m = l/(y/\n) is the 
bubble wall thickness. 

In two bubble collisions, the phase will at first -well 
before the bubble walls actually collide, at least in a zero 
temperature situation - interpolate between the values 
that it takes in each bubble by means of a phase wall 
situated at the midpoint in between the bubbles. This 
can be seen by taking the value of the modulus of the 
field p far away from the bubbles as that given by the 
ansatz 



^(bubblel + bubble!) = pe i9 ~ 

<&{bubblel) + <f>(bubble2) = Pl e lf>1 



Pie. 



Note that we will assume that (|7|) yields only a correct 
approximation for p, and not for 9. For p\,pi we will 
take an exponentially decaying ansatz in a 1 dimensional 
approximation 



Pi 
Pi 



Pt v e 

PtvC 



(x-\-vt — xq) / 5 m 
- (x— vt+x$) / Sjj 



(8) 



where v is the terminal velocity of the walls, Vt er above, 
and the bubble centers are situated at ±xq, with R < 
xq, R » I xo — R\ 3> 6 m (R being the bubbles radii). 
That is, we are basically asking that the distance from 
the bubble walls to the origin | Xq — R\ be several times S m 
(so that we fields decay exponentially there), and to have 
thin wall bubbles with R large enough so that the 1-d 
approximation is valid. Then (m) yields for the combined 
modulus 

p 2 ~ ^ e 2(„t-x )/5 m {2cosh(2x/(5 m ) + 2cos(6»i - 9 2 )} . 

(9) 

We can now use this expression in the equation of motion 
for the phase 9, to get 



2v 

Sm. 



2S- 1 sinh(2a;/£ n 



cosh(2a;/(5 m ) 



9' = 0, 



cos 



(10) 



where 6' = d9/dx,9= d9/dt. Note that for an initial 
phase difference 9\ — 9-2 — tt the modulus of the field 
(0) is zero at the midpoint between the bubbles. In this 
case then the denominator in the last term of ( |l0| ) goes 
to zero due to the fact that the phase has the shape 
of a step function as we go over the origin, switching 
discontinuously from 9\ to #2. On the other extreme, 
if the initial phase difference is zero there is of course 
no dynamics to it. We will then focus in an intermediate 
situation, and find solutions to ( |l0| ) for a phase difference 
of 7r/2. Using the ansatz 9 = T(t)X(x) 1 we can separate 
variables to obtain: 



T+ -=- T+ k 2 T = {), 
X" + A tanh (|^)X' + k 2 X = 0, 



(lla) 
(lib) 



k being the separation constant. The equation for T 
is obviously that of a damped oscillator. For k = 0, a 
T = const, solution exists consistent with our bound- 
ary condition that 9 goes to #1.2 as x — > =Foo at all 
times, while, for k 7^ all solutions will present unwanted 
temporal dependence temporal dependence of 9 even as 
x — > =poo. (Note that in case the lifetime of any transient 
configuration that might appear is only of the order S m /v, 
since the equation for T is that of a damped oscillator). 
We take k — then, whence an immediate integration 



of ( p.lb|) leads to the equation for X' (or equivalently 9' 
since T = const.) 



K 



cosh (2x/S m ) ' 



(12) 



with K being an integration constant with dimensions 
(length) -1 . This equation admits a family of general so- 
lutions given by the static sine-Gordon kink form 

9 = K5 m arctan(exp(2a;/<5 m )) + D = 

Carctan(exp(2ai/<5 m )) + D, (13) 

with K5 m = C, as can be easily checked by direct dif- 
ferentiation. (Note that this is not the original form 
in which these solutions appeared in the expression 
found there is equivalent to this standard sine-Gordon 
form though). The solution to the equation of motion for 
the phase (|10| ) interpolating from say 9% = at x — > — 00 
to #2 = at x — > +00 will then simply be 



2x 

1 (x) = arctan(exp( — )), 

Om 



(14) 



which clearly shows the structure of a phase wall placed 
at the origin and of width closely related to the bubble 
wall's width S m -although this relation will change as the 
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phase difference between the bubbles in (fO) change. It is 
also clear from this expression that this wall interpolates 
from 9\ to 9 2 following the shortest (geodesic) path in the 
vacuum manifold, since 9 monotonically increases from 
to tt/2. 

Since 9 admits sine-Gordon solutions though, it must 
be possible to write its equation of motion ( |l2| ) in a sine- 
Gordon manner. Let us use the form of 9(x) given by Jl3| ) 
and assume without loss of generality that the phase at 
x — > — oo is zero, so that D = 0. It is then not difficult 
to see that 




1 — cos 4 



Therefore, with V{9) 
sion 

d6 



)(x) 
C 



c 



1 



S m cosh(2x/i5 m )' 



(15) 



_ 



[1 — cos (4:9 /C)] , the expres- 



(16) 



is a sine-Gordon equation which has ( |l3| ) as its general so- 
lution. V has its minima located at 9 = mCir /2, m being 
an integer. Note though that (|l6|) is a static sine-Gordon 
equation, thus it only contains the kink and antikink sec- 
tors (apart from the vacuum of course) and its solutions 
only connect adjacent vacua. 

Going back to our particular problem, where A9 = 
9 2 — 9 1 = n/2, we can now use the fact that the physical 
phase 9 is defined only up to modulo 2tt to choose 9 2 
to be 9 2 = tt/2 + 2nn, with n = ... - 2,-1,0, 1,2, .... 
being an integer. It immediately follows that (|l4|) is not 
the only solution consistent with our physical boundary 
conditions. In fact there is an infinite family of solutions 
to ( [l2] ) having A9 = tt/2 + 2irn which will satisfy our 
requirements, of which ([u]) is only the n = case. An 
interesting case corresponds to n = — 1. This phase wall 
forces 9 to monotonically decrease from 9\ = at x — > 
—00 to 9 2 — —3tt/2 at x — > -f-00, and therefore follows 
the shortest non-geodesic path in the vacuum manifold 
between 9\ and 9 2 . Its form is 



2x 

9 (x) = — 3 arctan(exp( — )). 



(17) 



Lower negative values of n will correspond to phase walls 
that completely wind around the vacuum manifold — \n\ + 
1 times in the negative 9 direction before going from 
to — 3tt/2, whereas positive values of n will lead to phase 
walls carrying a (positive) winding number n. (Hereafter 
in the paper, when using the word "winding" to denote 
the phase change across a phase wall it should always be 
understood that there is no phase circulation around any 
closed path in physical space, as opposed to the times 
in which "winding" is used to refer to the circulation of 
the phase around a vortex or antivortex). The expression 
giving 9(x) for a general value of n will then simply be 



2x 

9 (x) = [1 + 4ri] arctan(exp( — )). 



(18) 



III. PHASE WALLS IN 2 SPATIAL DIMENSIONS 
WITHOUT DISSIPATION 

We will now extend the analysis carried out previously 
to the case in which there is no dissipative coupling of 
the scalar field to the rest of the matter present in the 
model. Also, and since we will be needing it in order 
to be accurate in the numerical simulations that we will 
perform, we will introduce 2 spatial dimensions in the 
problem. 

The problem of two bubbles approaching each other 
and colliding in 2+1 dimensions is invariant under the 
SO(l,l) Lorentz group. If we take the bubbles to nu- 
cleate initially at (0,±xo,0) with coordinates (t,x,y) and 
signature (—,+,+), then the symmetry of the problem 
dictates that we must have p = p(x,t) and 9 = 9(x,t), 
where r 2 = y 2 — t 2 . Because of the particular expression 
that we will use as the approximate value of the modu- 
lus of the field though, a different set of coordinates turn 
out to be more convenient. We define new coordinates 
u, v, a 



Note that 



x = xq cosh u cos v , 

y = xq sinh u sin v cosh a, 

t = xosinhusinwsinha. 



t = xq sinh u sin v, 



(19a) 
(19b) 
(19c) 



(20) 



thus to preserve Lorentz invariance in these coordinates 
is simply to say that p and 9 can not depend on a. The 
range of variation of u,v,a will be 

u > 0; -tt/2 < v < tt/2; a > 0, (21) 

and in terms of these coordinates the equation of motion 
for 9 will read (dropping already any dependence on a) 

d U u9 + d vv 9 + coth u- d u 9 + cot v ■ d v 9 

+- (d uP d u 6 + d vP d v 9) = 0. (22) 
P 

Now we have to specify an ansatz for p. As in the case 
before, we will take the modulus of the field between the 
bubbles to be given by the modulus of the sum of the 
two individual bubble fields. Thus, 

|$| 2 =p 2 = \ Pl e^ +p 2 e^\ 2 =p 2 +p 2 + 2p lP2 cos(A9). 

(23) 

Let us now define the invariant distance of any point to 
the bubble nucleation events, xi,2 



Xl,2 



x ± xo) 



V 



t 2 . 



(24) 



where the plus sign holds for ^and the minus for \2- 
Looking back at (|]) it is obvious that outside the bubbles 
and away of the their walls we should have 
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Pi, 2 — Ptv& 



-(Xl,2-Ro)/8 n 



(25) 



Taking this into (|23|), setting already A9 = n/2, and 
writing the result in terms of u and v leads to an expres- 
sion for p 



P 



•r, 2 "» (^cosu) . (2(i! 



studied in the previous section, n = will correspond to 
the geodesic phase wall solution, while higher (or lower) 
values of n lead to phase walls that wind around the 
vacuum manifold of the model n (or \n\ — 1) times. As 
before, n = — 1 corresponds to the shortest non-geodesic 
path. We have thus shown that the existence of phase 
wall solutions is a general feature of the model. 



as can be checked by trivial, if lengthy, manipulations. 
Then plugging this form for p into (|22|) results in 



coth u — - sinh u 



d U u9 + d vv 9 + 

2x . / 2x n 

cot v — 7 — sin v tanh — — cos v 

Om. \ Om 



d v 9 = 



(27) 



In spite of the appearance of this equation, it is actually 
not too difficult to find the type of solutions that we are 
looking for, at least within some approximation. Note 
that we expect the phase walls to have their central axis 
located at x = 0, the line equidistant to the bubble cen- 
ters, which in u,v coordinates corresponds to v = tt/2. 
If we expand the second bracket in powers of v around 
vq = 7r/2 it is easy to see that neglecting the cotu term 
when compared to the other term in the bracket amounts 
to neglecting terms of order 1 or smaller when compared 
to terms of order (2xo/S m ) 2 or higher in the coefficients 
of the series. Since within the approximation that we are 
working with xo/S m 3> 1, we conclude that to an accu- 
racy 1/ (2xo /o~m) 2 we can safely neglect the first term in 
that bracket. Assuming as ansatz that the solutions that 
we are looking for depend only in v leads then to 

d vv 9 j—— sin v tanh ( -—— - cos i> ] d v 9 = (28) 

Om V °m J 

which indeed has phase wall solutions of the form 



= [An + 1] arctan(exp(—^ cosv)), 

Om 



or, in terms of xi 2 



9 = [An + 1] arctan(exp( 



Xi - X.2 



))• 



(29) 



(30) 



Expression ( 30|) is clearly again a sine-Gordon kink, dif- 
fering from ( 18[) only in its argument. Note however that 
the fact that (|30|) depends on \i~ Xi implies that it does 
not represent a planar wall. This is depicted in figure 
[I], where we show constant phase lines for the geodesic 
phase wall interpolating from to 7r/2. Thus, lines of 
constant phase will have a curvature which will depend 
on the distance to the bubble centers (in the figure the 
bubble centers are situated at xq — ±10). If we expand 
Xi,2 m terms of ~X*/xq though, where ~x* = (t,x,y), ( |30| ) 
reduces to the 1-dimensional expression ( |l8| ) above for as 
long as (x/xo) 2 ~ (y/xa) 2 ~ (t/xo) 2 -C 1. As in the case 



IV. PHASE WALL STABILITY ANALYSIS 

Since an analytical approach proved to be too involved, 
we resorted to numerical simulations in order to investi- 
gate the problem of the stability of non-geodesic phase 
walls -the geodesic one being obviously stable. We used 
a straight forward leapfrog method with second order 
accuracy to evolve in time different initial configurations 
corresponding to phase walls of various winding numbers. 
In order to better understand the process of phase wall 
decay we first studied the evolution of the 1-dimensional 
solutions, then carried on the study to two dimensional 
configurations. 

We begin in figure || by plotting, using polar coordi- 
nates in field space, the fate of the shortest nongeodesic 
wall through a series of snapshots. In||(a) we can see the 
initial configuration of the phase at t = 0, going from 
9 = 0to9 = — 3n/2. Subsequently we see how in the very 
core of the wall, the field starts decreasing its modulus, 
eventually crossing the origin (p = 0), and finally settling 
down in the geodesic configuration, in which 9 goes from 
9 = to 9 = it/2, by t = 2 in units of the bubble wall 
thickness S m . In figure p| we have depicted the unwind- 
ing of an = —2 wall. In0(a) we can see how the initial 
configuration has, as expected, an extra 27r winding over 
the n = — 1 wall previously seen. Note how in this case 
the decay proceeds by unwinding first this extra 2-7r loop, 
after which we are left with a configuration equivalent 
to that of a n = —1 wall which in turn decays as de- 
scribed above. The unwinding of the extra 2-7T is shown 
in figure ^. In it we can see how the loop is progressively 
displaced and shrunk, until the field crosses the origin 
at the midpoint between the bubble centers. After that, 
the field reeqnarrays itself. The gradient energy acts as a 
tension on the phase loop, which eventually evolves into 
a cusp. This cusp then divides into two sharp kinks that 
move away from each other. As we can see then the main 
features that seem to emerge of these results are that a 
phase wall of winding n will cross the zero field value \n\ 
times in its decay, decreasing n in steps of one, and that 
the time scale of the decay is of the order of 5 m , becoming 
progressively shorter as the winding increases. 

Using these results it is not too difficult to convince 
oneself of the fact that in two spatial dimensions we 
should in general expect the decay process to produce 
pairs of vortices and antivortices. Let us think for in- 
stance of the symmetry axis of a n = —1 phase wall, 
that is, the line of points that are equidistant to the bub- 
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ble centers. (For a 1-dimensional wall, the symmetry 
axis is of course reduced to a single point: the central 
point equidistant to the bubble centers). In two space 
dimensions we may think of the n = — 1 wall decay as 
being the result of all the points in this axis crossing 
the zero field value -each following a similar process as 
the 1-dimensional case- and reeqnarraying to follow the 
geodesic path. The points along the wall's symmetry 
axis are not all equivalent to one another though, since 
the distance to the bubble centers changes as we move 
along the axis. Therefore the field does not cross the 
zero value everywhere along the axis simultaneously, but 
rather some points will go through this value before oth- 
ers. In particular, given the symmetry of the problem 
it seems clear that the central point of the axis, with 
the smallest distance to the bubble centers, will cross the 
zero first. Once the field has switched to the geodesic 
path along the line that joins the bubble centers, the two 
points immediately above and below the central point 
will cross the zero themselves, then the ones respectively 
above and below these two, and so on. Therefore, as the 
phase shifts towards the geodesic path we always have a 
pair of points on the wall axis at which the field takes the 
false vacuum value, both of them quickly moving up and 
down respectively. Since the direction of the field in the 
region that has already switched to the geodesic path is 
opposite to its direction immediately above and below it, 
we end up having a net circulation of the phase around 
these two points. We have thus produced a vortex- 
antivortex pair along the phase wall axis, the two defects 
moving away from each other. 

In figure [| we can see how this process takes place. 
In it we can see the contour lines of the modulus of the 
field, while the phase is represented by the angle that the 
arrows form with the x axis. The length of the arrows 
is proportional to the modulus of the field at the point 
at which it is drawn, except in the central region with 
p < 0.001 where they have been magnified to a com- 
mon length so that we can appreciate the decay process. 
Throughout the simulations the true vacuum value of the 
field is ptv ~ 1. In figure [^(a) we have the initial n = — 1 
configuration, if we follow the change of the phase hori- 
zontally along the central y = line for instance we can 
see how the phase steadily decreases from to — 3ir/2. In 
figure ||(b) we see how the field starts decreasing its mod- 
ulus, especially in the central region between the two bub- 
bles. Eventually, as described above, the central point 
will cross the p = value first. After this has occurred 
the points above and below it will follow suit, then the 
ones respectively above and below these two, and so on. 
Due to the opposing phase directions above and below 
the p = points though a vortex-antivortex pair will 
form, as can be seen in||(c). As successive points along 
the axis cross the zero, the defects move further and fur- 
ther away from each other, as shown in ||(d). Note how 
on |^(d) the phase still has some reeqnarrayment left to 
do in the region that has crossed the zero field value. It is 
below zero right before the wall axis, and over 7r/2 right 



after it as we move from left to right. This situation is 
analogous to what we see in figure g(b) and (c) in the 1- 
dimensional case. Eventually, the geodesic configuration 
is attained. 

In figure ^ we can see the same process taking place 
for the n = — 2 phase wall. As was the case with the 
1-dimensional solution, this configuration first decays by 
unwinding its extra 2n to become an— —1 wall. In or- 
der to do this however it must cross the p = value as we 
saw in 1 dimensions. The 2tt unwinding will thus make 
the field cross its zero value, first at the origin, then gen- 
erate two zeroes of the field that move away from each 
other along the wall axis following a process analogous 
to the one just seen. As before thus, a pair of defects is 
created in this decay. Although not shown in the figures, 
the n = — 1 wall will then in turn decay to the geodesic 
configuration generating a second pair of defects. There- 
fore, the final unwinding of the n = —2 wall produces 
two pairs of defects. The rule that immediately seems to 
suggest itself is that a phase wall of winding n will pro- 
duce precisely |n| pairs of defects in its decay. The simple 
reason behind it would be that for such a wall the field 
needs to cross \n\ times the p = value in order to com- 
plete its decay, and that each time that the field crosses 
the zero value at the central point a vortex-antivortex 
pair is produced. One should be careful though. Note 
that the decay of the n — — 2 wall to the n = — 1 config- 
uration sends and vortex into the lower half axis of the 
wall, whereas the decay of the n = — 1 wall sends an an- 
tivortex into the same region. Since global vortices and 
antivortices attract each other one could worry that these 
pairs could tend to annihilate, especially in the case of 
high n walls where the defects will be produced within a 
short distance of each other. It would seem though that 
such an annihilation can not actually happen, since each 
defect carries away as phase circulation around it one of 
the winding factors present in the wall. (I.e., think for 
instance of the n = —2 case. The successive vortex and 
antivortex produced in the decay can not annihilate each 
other because, if they did, the phase would have to go 
from the region in which the wall has not decayed yet 
into the region in which it has, having no zeros around 
which it can wind in between). It is unclear at the mo- 
ment whether this same argument holds for higher n walls 
however. More simulations should perhaps be carried to 
test this case. 



V. NON GEODESIC PHASE WALLS AND 
THREE BUBBLE COLLISIONS: TOPOLOGICAL 
DEFECT NUCLEATION 

If non geodesic phase walls can indeed be formed in two 
bubble collisions and they subsequently decay by nucle- 
ating vortices, then the next obvious question is: what 
will be the impact of these processes on the probabil- 
ity of topological defect becoming trapped in a region of 
true vacuum when three bubbles collide?. The standard 
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picture of defect nucleation described in the introduction 
roughly yields a prediction of 1 defect per every three 
bubbles (times a proportionality constant to account for 
the probability of the three bubbles having the appropri- 
ate phase). Let us think of three colliding bubbles with 
phases such that they would nucleate a defect following 
the geodesic rule: is the defect still formed when the 
phase chooses a non-geodesic path between two of the 
bubbles?. Alternatively, if the initial phases would not 
yield a vortex following the geodesic rule, will a vortex be 
nucleated if a non-geodesic wall is formed between any 
two bubbles?. Clearly, these questions can potentially 
alter the standard prediction for vortex nucleation prob- 
ability. We have therefore simulated these two scenarios 
to find out their answers. 

In figure [7] we can see the results of the first case. In it 
we have three bubbles with phases 7r/6, 57r/6 and 3ir/2 
that in principle would nucleate a vortex following the 
geodesic rule. The two upper bubbles however have col- 
lided slightly before, hitting the third one a little later, 
and the phase has wound up between them following a 
n = — 1 phase wall. As we can see in figure ^(a) though, 
the field does not reach its zero value anywhere in the 
region between the three bubble centers, as it would be 
the case if the phase were to follow the geodesic path 
everywhere. This is the key of the question here: if the 
geodesic rule were to be holding, the vortex predicted by 
it would appear in the region between the three bubbles 
as soon as they make contact with each other, in this 
case however the vortex is not there to begin with. As 
we have seen though, the phase wall has to decay pro- 
ducing a vortex-antivortex pair. In the following figures 
we can then see how the vortex produced in the decay 
ends up resting between the three bubbles, while the an- 
tivortex escapes "upwards" and away from the bubbles, 
into a region that gets exponentially closer and closer to 
the false vacuum. Thus, in terms of stable vortices -that 
is, vortices surrounded by regions of true vacuum- the 
net result is exactly the same as the one predicted by 
the geodesic rule. There is however, one extra antivortex 
that has been expelled into the near-false vacuum region. 

The second case is depicted in figure ||. Clearly the 
phases inside the bubbles would not produce a vortex 
should the geodesic rule hold. As before, a n = — 1 phase 
wall has formed between the upper two bubbles, and in 
this case the field already has a zero in the region between 
the bubble centers as we can see in figure ||(a). When 
the phase wall unwinds, the resulting vortex annihilates 
with the existing antivortex, therefore yielding again the 
same result as the geodesic rule would predict: no stable 
defects are formed. Again, the other antivortex resulting 
from the phase wall decay escapes away from the bubbles. 

Although we have not carried out simulations with 
phase walls of higher winding, it is easy to convince one- 
self that the results will not change for higher |n| configu- 
rations. For instance, to reproduce with a higher winding 
an initial situation analogue to the one in figure ^(a) we 
would need a phase wall with n = —3. However, as this 



wall decays first to n — —2 and then to n — — 1 it nucle- 
ates into the region between the bubbles first a vortex, 
then an antivortex. Obviously this pair annihilates as it 
gets trapped between the bubbles, thus leaving us with 
the same initial state as the one we have studied. The 
only difference would lie in the number of defects that 
are expelled away from the bubbles. 

VI. A MODEL FOR HOW THE GEODESIC RULE 
WORKS IN PHASE TRANSITIONS WITH 
GLOBAL U(l) SYMMETRY 

Since the non-geodesic phase walls and their decay de- 
scribed in the sections above were found to form in the 
region between the bubbles where the modulus has an ex- 
ponentially low value, the preceding considerations would 
seem to be only an academic exercise. In this section we 
will however attempt to give a physical picture of how 
the geodesic rule could work in models with global sym- 
metry. 

We will start by assuming that the physical set up in 
which the process takes place is that of a first order phase 
transition at finite temperature. The oscillations in the 
value of the field around its false vacuum value due to the 
finite temperature effects will yield a background thermal 
field. These field fluctuations will each have a phase, and 
we will denote the mean value of their modulus over space 
and time as < p t h >• They will also have an average size 
(radius) and lifetime, < £ > and < r >. We will also 
assume that < p t h > is still well below the true vacuum 
value of the field p tv so that the thermal fluctuations do 
not strongly affect the field inside each bubble and the 
values of the phase of each of them, #iand 82, remain well 
defined and largely unaffected by the fluctuations. It is 
against this thermal background that the bubble fields, 
0iand 02, will propagate. The p\^ will obviously have 
the shape of the bubble profile, but any trace of them 
will be lost soon after their value becomes lower than 
< pth > since beyond the bubble walls they decay ex- 
ponentially with a typical decay distance S m . Therefore 
we will assume that beyond a few d m from the point at 
which pi.2 ~< Pth > the thermal background field will 
effectively erase any presence of the bubble fields. Figure 
^ exemplifies this picture. The two bubbles will first col- 
lide along the line that joins their centers. Roughly, this 
will happen at about the time when their combined mod- 
ulus at the point equidistant to the centers pb -given by 
the ansatz (^) evaluated at x — 0- reaches a value equal 
to the thermal fluctuation average pt ~< p t h >, and the 
whole collision process will take place on a region that 
we will characterize by a length d. Although d will most 
likely be of the order of a few S m , we will leave it largely 
unspecified for the time being. As we will see this will 
not have a big impact in our conclusions as far as this 
simple argument is concerned. When the collision takes 
place the phase has to decide which path will it follow 
in the vacuum manifold in order to interpolate from 9\ 
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to 02- Obviously the value of the phase of the thermal 
fluctuations in the collision region will have an important 
influence on what configuration is finally chosen, since by 
continuity the phase will be forced to interpolate between 
6\ and 82 going through the value that it takes at the col- 
lision region. After that, pb will rise above < p t h > expo- 
nentially fast and as far the thermal effects are concerned 
the configuration chosen by the phase at that point will 
remain frozen thereafter. We will work in three spatial 
dimensions, and as elsewhere in the paper we will take 
#i=0 and 62 = 7r/2. This case is of particular interest 
because since A6 will always lie randomly between and 
7r (modulo 2tt) , the average phase difference between any 
two given bubbles will precisely be AO = it/2. We will 
examine two limiting cases. 

As our first case let us examine a slow collision sce- 
nario in which the time scale of the collision t co i (i.e., the 
time needed by the bubbles to expand their radii by an 
amount (Q) is much larger than the average lifetime of the 
fluctuations < r >. If t co \ ^>< r >, fluctuations appear- 
ing between the bubbles when they are at a distance d or 
less away from each other will be energetically favored to 
have a phase that falls within the geodesic path from 9\ 
to 62- Under these conditions, we will suppose that the 
probability p„ of choosing any given configuration with 
winding n is weighted by a Boltzmann type of distribu- 
tion, p n ~ e ~P E ™ where as usual (3 = 1/T (in natural 
units, k = 1), and E n is the energy corresponding to the 
given configuration. Thus, the relative suppression fac- 
tor for non-geodesic walls of winding n with respect to 
the geodesic configuration will be given by e~ l3AE \ where 
AE = E n — Eo. Since both solutions have the same back- 
ground modulus for the field though, the gradient of the 
modulus does not contribute to this energy difference, 
and neither docs the potential energy since the potential 
depends only on the modulus. If we then assume that 
the collision region is small enough so that the three di- 
mensional energy density can be approximated by its one 
dimensional limit, the total energy difference will simply 
be given by the one dimensional integral times the cross 
sectional area of the collision region, d 2 . As was seen at 
the end of section II this approximation will be valid for 
as long as {d/xaf -C 1, with xq being the distance to the 
bubble centers, which will be justified in our case. Thus, 
as an order of magnitude estimate we will have 

/+00 
pI (del - del) dx ~ 
-00 

/ + OO 
(del - del) dx, (31) 
-00 



*For bubbles expanding at the speed of light t co i ~ d, note 
however that if the motion of the bubble walls is damped 
by interactions with the external plasma (as the case seen in 
section I) t co i can be significantly larger than d. 



since the phase gradients are exponentially close to zero 
outside of the phase wall, where at the moment in which 
the collision takes place pb ^< pth >■ We can now use 
the well known results for the sine-Gordon kinks to get 

AE ~ <Pth >2 d 2 [{An + If - 1] . (32) 

This is the key expression that will tell us how suppressed 
the non-geodesic configurations will be. We can see how 
phase walls with high winding get quickly suppressed, as 
one would have expected. As expected too, formation 
of non-geodesic walls is also strongly suppressed for suf- 
ficiently high values of < p t h >■ This agrees with the 
standard argument which states that at least for transi- 
tions with global symmetry the geodesic rule must hold 
since it yields the least energetic configurations. Note 
however how, for sufficiently low values of < pth > this 
argument doesn't hold any more and production of non- 
geodesic walls could take place, in particular those with 
low winding. The underlying physical reason is simply 
that the energy density of the phase wall is given by 
p 2 d0 2 , therefore, if the wall is placed in a region were p 2 
is small the energy difference between the different con- 
figurations will also be small. Once the phase has chosen 
an interpolating path, the modulus of the bubble field 
will quickly rise against the thermal background, and the 
influence of the thermal fluctuations will become negli- 
gible as stated above. The non-geodesic walls will sub- 
sequently decay, the pairs of defects generated by their 
decay being expelled into the thermal background. It 
would seem though that in this case we can expect to 
form non-geodesic phase walls only for small values of 
< Pth >■ 

As our second limit let us examine the fast collision 
scenario and suppose that t co i <C< r >. In this case 
the last fluctuation on the scale of the collision region 
appeared well before the bubbles were close to it. Let us 
then make the naive assumption that, largely unaffected 
by the bubbles, fluctuations in the collision region pick 
up phase values e p at random from to 2-7T at the time 
at which they nucleate, with all values having the same 
probability. Since by definition the geodesic path from 0\ 
to 62 is the shortest one on the vacuum manifold, the field 
fluctuations in the collision region are more likely to have 
a phase that falls out of the geodesic path than inside of 
it. If, whatever value chosen for p , the paths travelled 
by the phase to interpolate from 0\ to e p and then from 
e p to 6>2 are geodesic, then it is easy to see that any value 
of Op between tt and 37r/2 will lead to the formation of 
a n = — 1 phase wall. That is, the probability p(ng) of 
forming a non-geodesic phase wall in this case will be 
simply 

Ping) = ^ = \. 03) 

Note that this is an estimate from below, since phase 
walls with winding other than —1 are not being consid- 
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ered at all in this simple picture. (They would corre- 
spond to non-geodesic interpolations between the bubble 
phases and 6 P ). Note also how in this limit the fact that 
for any two bubbles A9 = ir/2, as pointed above, also 
implies that the average probability of forming a non- 
geodesic wall between any two bubbles will be 0.25, since 
the fact that < AO < it with all values having equal 
probability automatically implies that < p(ng) < 0.5 
(we exclude the case AO — tt since in that case any of 
the two paths is geodesic), and that p(ng) — 0.25. As 
before however, any nongeodesic phase wall will quickly 
decay into the geodesic configuration, expelling vortices 
and antivortices into the background. 

The simple argument we have just presented would 
seem to indicate then that in the case of slow collisions 
the geodesic rule is established mostly by the thermal 
suppression of non-geodesic configurations, while for fast 
collisions the rule holds because non-geodesic walls are 
unstable and decay into the geodesic one. Since in both 
limits the geodesic rule holds, one would also expect this 
to be the case in the (more physical) intermediate regime 
in which t co i << r >. 

VII. CONCLUSIONS 

Although the geodesic rule is known to hold for first 
order phase transitions that break a global U(l) symme- 
try, a dynamical account as to how the field settles into 
the geodesic configuration was so far missing in the liter- 
ature. In this paper we have attempted to give precisely 
such an account, in the hope that some of the findings 
and the methods here developed could shed some light 
on the more difficult questions concerning the validity 
of the rule in the Abelian-Higgs model. Using a fam- 
ily of analytical solutions that represent the phase as it 
interpolates between its values inside two different bub- 
bles, i.e., phase walls, we have found both the solution 
corresponding to the geodesic path as well as infinitely 
many more which wind around the vacuum manifold n 
times before reaching the desired value. After finding 
the expression for the phase walls in 2+1 dimensions, we 
turned our attention to the question of the stability of the 
non-geodesic configurations. Our analysis shows that, as 
one would expect, non-geodesic phase walls are unstable, 
decaying into the geodesic configuration in times of the 
order of S m . 

The decay of the non-geodesic states takes place via the 
nucleation of vortex-antivortex pairs that quickly move 
away from each other, and the simple rule that seems 
to be the result from our study is that a phase wall of 
winding n will produce precisely |n| pairs of defects in its 
decay. Although this production of defect pairs does not 
alter the geodesic rule standard prediction of stable de- 



fect[] nucleation in three bubbles collisions, it is perhaps 
worth pointing out that for first order phase transitions 
with supercooling (where < p t h > will be low and < r > 
large) a significant number of bubble collisions may po- 
tentially result in the production of non-geodesic phase 
walls and their subsequent decay. Thus, it would seem 
that for these transitions one could think of the thermal 
background as of a region in which a gas of defects and 
antidefects exist. Since the average modulus < p t h > 
of this thermal field will be low in supercooled transi- 
tions, the origin of these defects could in general be due 
to both thermal fluctuations of the field and to their pro- 
duction in phase walls decay. For as long as the second 
mechanism is not dominant, one could think of this gas 
as being roughly in equilibrium, since, in a large volume 
thermal fluctuations would be creating and annihilating 
defects and anti-defects with roughly the same probabil- 
ity. If production via phase wall decay becomes dominant 
though, the density of defects in the thermal background 
will tend to grow as the transition proceeds. A number of 
questions are raised by this scenario. For instance, what 
will be the final fate of all these defects as the transition 
reaches its end?. Will they all annihilate one another 
or will a significant fraction survive the transition?. Al- 
though in a large volume the average number of vortices 
will be equal to that of antivortices, in a single smaller 
volume vortices and antivortices need not exactly can- 
cel each other out. What will happen if a three bubble 
collision traps a net winding of the phase in the basin 
between the bubbles?. Also, could the presence of this 
defect gas in the thermal background have an impact on 
models in which baryogenesis is induced by topological 
defects?. 

As far as the connection with the Abelian-Higgs model 
is concerned, some of the results that we have found could 
be of interest. The motivation for our study lies partly 
in the hope that one could export to the gauge case the 
same type of techniques that we have developed to study 
the global model. This could be done by defining some 
gauge invariant form of the phase difference (i.e., a la 
Kibble- Vilcnkin for instance ||, where A9 = J dx k DkO, 
with D^O — d^O + eA^, and e, being the gauge cou- 
pling constant and gauge field respectively) and studying 
its equations of motion. It is also worth pointing out the 
close resemblance that some of our results bear to those 
of 0. In this work, the authors found out that both in 
the case of the Abelian-Higgs model and the sin 0\y — 
limit of the Weinberg-Salam model there may be enough 
energy released from the collision of the bubble walls as to 
drive the modulus of the Higgs away from its vacuum ex- 
pectation value and make it oscillate around zero. Extra 
regions of symmetry restoration are thus produced, and 
the field oscillations as they decay lead to the formation 



i.e., defects surrounded by regions of true vacuum 
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of multiple windings. In their simulations, they found 
that the validity of the geodesic rule seemed to depend 
not only on the gauge coupling, but also in the initial 
bubble separation and the phase difference between the 
bubbles. Although in our model the energy needed to 
take the phase wall away from the geodesic configura- 
tion is provided by the background thermal fluctuations, 
rather than wall collision, note that our two dimensional 
phase walls also show a dependence both on the distance 
between the bubble centers and on the initial phase dif- 
ference. (Although only for A9 — tt/2 do we have an 
analytical solution, solutions for other values of A8 are 
in principle obtainable numerically). Furthermore, high 
winding phase walls decay via field fluctuations around 
the false vacuum value, nucleating multiple defects as 
they decay. 
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FIG. 1. Lines of constant 9 for a geodesic phase wall going 
from 9 = to 8 = 7r/2. The lines are labelled in units of A6. 
Thus, a label of 0.5 implies a value of the phase of 7r/4, and so 
on. The bubble centers are placed at xq — ±10, and we have 
taken 8 m — 1. A lower value of |xo| would make the contour 
lines show more curvature for the same axis scale, whereas for 
larger values the wall would look flatter on this scale. 

FIG. 2. The decay of a n — —1 phase wall is shown in 
polar coordinates in field space. Distance to te origin gives 
the modulus of the field, while the angle corresponds to the 
field phase. In (a) the initial configuration at t = is shown. 

(b) corresponds to t = 0.4, (c) to t = 0.6, and finally (d), 
in which we can see that the geodesic wall is the final state, 
corresponds to t = 2. Time is measured in units of the bubble 
wall thickness 5 m - 

FIG. 3. We see here the decay of the n = — 2 phase wall 
into a configuration equivalent to a n — —1 wall. As before, 
the initial t = configuration is shown in (a). Subsequent 
images correspond to: (b) t = 0.1, (c) t = 0.15, and (d) 
t = 0.2 . 

FIG. 4. The shrinkage of the extra 2tt loop is shown here in 
greater detail. Images correspond to: (a) t = 0, (b) t = 0.11, 

(c) t = 0.13 by which time the field has already crossed the ori- 
gin and we no longer have a 2-k phase change around the loop, 

(d) t = 0.14 where we see how the loop has been stretched 
into a cusp, (e) t — 0.16 where the cusp has divided into two 
sharp kinks that move away from each other, and finally (f) 
t = 0.18 . 

FIG. 5. A n — — 1 two dimensional phase wall decays into 
the geodesic configuration. Solid lines indicate contour lines 
of the field modulus, while the phase is given by the angle 
between the arrows and the x axis. The length of the arrows 
is proportional to the value of the mudulus at the point they 
are drawn, except in the central region occupied by the phase 
wall, wher they have been magnified to a common length. The 
initial configuration at t = is shown in (a). In (b) t = 0.3, 
and we see the modulus of the field quickly decreasing in the 
central region. In (c) t = 0.304, the field has crossed the 
origin in field space at the central point, and its modulus is 
now growing there. As described in the text this produces 
a vortex -antivortex pair, we can see the vortex somewhat 
above y = 1, and the antivortex at y = — 1. In (d) t = 0.308, 
the defects keep moving away from each other as the field 
increases its modulus in the central region. 

FIG. 6. A 2-dimensional n — — 2 phase wall decays into 
the n = — lconfiguration. The initial coniguration is shown 
in (a) at t = 0. To see the extra 2n factor in the phase 
follow its change along a horizontal line, the arrows turning 
always in a clockwise direction. In (b) t = 0.108, the field 
is decreasing its value in the central region, in (c) t = 0.109 
and the vortex (around y = — 1) and antivortex (y — 1) have 
already appeared. Note how between the defects the phase 
is now increasing across the wall from almost — tt to almost 
—tt/2, subsequently decreasing again until it reaches — 37r/2. 
In (d), t — 0.111, the defects keep moving away from each 
other, and the field increases its modulus in the central region. 
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FIG. 7. a n — — 1 phase wallbetween two bubbles with 
phases 5n/6 (left) and 7r/6 (right) decays in the presence of 
a third bubble with phase 3tt/2. In (a) t = 0, (b) t = 0.34, 
(c) t = 0.384, and (d) t = 0.8. The decay sends a vortex into 
the region between the three bubbles that gets trapped there. 
The (upper) antivortex is expelled away from the bubbles. 



FIG. 8. A n = — 1 wall decays in the presence of a third 
bubbe. The three bubbles would not form a defect accord- 
ing to the geodesic rule. The net result of the process agrees 
with that predicition. In (a) t = 0, we can already see how, 
as a consequence of the presence of the n = — 1 wall we al- 
ready have an antivortex in the region between the bubble 
centers. The wall subsequently decays, in (b) t = 0.34, and 
(c) t = 0.372, sending a vortex into the interior region. This 
vortex and the antivortex annihilate each other, leaving in (d) 
t = 0.8 the final configuration predicted by the geodesic rule. 



FIG. 9. Two bubbles bubbles about to collide expand 
against a fluctuating thermal field background. Since the bub- 
ble fields decay exponentially fast beyond the walls, any trace 
of them is quickly lost once they get below the thermal field 
level. The thermal fluctuations are presumed not to be strong 
enough to significantly alter the bubble fields inside the walls 
however. 
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